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The AH(M, AM) method and its ability to determine intrinsic switching field distributions of 
perpendicular recording media are numerically studied. It is found that the presence of dipolar 
interactions in the range of typical recording media substantially enhances the reliability of the 
AH(M, AM) method. In addition, a strong correlation is observed between the precision of this 
method and a self-consistency-check of the data sets, which is based upon a simple redundancy 
measure. This suggests that the latter can be utilized as an efficient criterion to decide if a complete 
data analysis is warranted or not. 

PACS numbers: 75.50.Ss, 75.60.Ej, 75.75.+a 
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For magnetic recording media used in state-of-the-art 
storage applications such as hard disks drives, the in- 
trinsic switching- field distribution D(Hq) of the media 
grains is one of the most crucial properties defining the 
recording quality [![. In general, each grain is charac- 
terized by an intrinsic switching field H$, which is a lo- 
cal material property. Because grains interact with each 
other by means of exchange and dipolar interactions, 
D(Hs) is not easily accessible in macroscopic measure- 
ments. This is especially true for perpendicular record- 
ing media due to the strength of the interactions, which 
is much larger than in the previously used longitudinal 
recording media Therefore, many attempts to de- 
termine D(Hs) from macroscopic magnetization rever- 
sal type measurements have had varying success with 
none of the presently available methodologies being com- 
pletely satisfactory for all levels of intergranular interac- 

tionsiinaiai. 

The recently developed AH(M, AM) method has been 
used in analyzing and quantifying progress in perpendic- 
ular recording media fabrication @j U 9 • This method 
measures the field difference AH at constant magnetiza- 
tion M between the major hysteresis loop and a number 
of recoil curves, which each start at a certain distance 
AM away from saturation. Within the mean-field ap- 
proximation, the functional dependency of AH can be 
written as AH(M, AM) = - M)/2) - - 

M — AM)/ 2) where J -1 is the inverse of the integral 
I{x) = f_ D(Hs)dHs- Within the framework of this 
method, AH is independent from the grain interactions, 
which allows for a direct experimental access to determin- 
ing D(Hs). For certain parameterized distribution func- 
tions, one can derive analytic expressions for AH. For 
example, for a Gaussian distribution of width a, one finds 
AH G (M,AM) = y/2a[ed~ l (M+AM)-erT 1 (M)]. De- 
tails of this method and the analysis formalism have been 
described previously @, @, H, [t| . 

The AH(M, AM) method was demonstrated to have 
several advantages over comparable methods. First, it 



allows for the determination of the entire D(Hg) distri- 
bution and its functional form and not just a single char- 
acteristic parameter 0, H[ . Second, it has a well-defined 
reliability range and it allows for oversampling, which 
makes self-consistency checks feasible @. Third, its fail- 
ure mode was found to show universal behavior, inde- 
pendent from the detailed lattice structure in numerical 
simulations [Toj |. 

Despite these advantages, it is still unknown whether 
the reliability range of the AH(M, AM) method with 
respect to exchange interactions might be affected by the 
simultaneous presence of dipolar interactions. Given the 
fact that this simultaneous presence of both interactions 
is the realistic case for actual recording media, it is an 
important issue, which we have studied in this letter. 

For our numerical studies, we model each media grain 
as a symmetric hysteron, which generates a rectangu- 
lar hysteresis loop in an applied field H Q. The half 
width of the hysteresis loop is just the intrinsic switch- 
ing field Hs of this hysteron. We further assume that 
the magnetization of each hysteron exhibits values of 
±1 only and orients exactly along the applied field H. 
The ferromagnetic layer system is then represented by 
a square or triangular lattice of symmetric hysterons 
with periodic boundary conditions. The model Hamil- 
tonian can then be written as Ji = — J C xJ2<i j> $iSj + 

^dpE^i ff 1 - E*(H + B^{Si)H Si )Si. Here, the first 
term represents that hysterons interact ferromagnetically 
with their nearest neighbors by means of exchange in- 
teractions of strength J ox . The second term represents 
that hysterons exhibit a distance-dependent dipolar in- 
teraction of strength J^p with all other hysterons. The 
third term accounts for the effect of the external field 
and the intrinsic switching field. This model is referred 
to as the interacting random hysteron model. The algo- 
rithm for the simulation of the major hysteresis loops and 
the corresponding recoil curves of this model is described 
in Q . For the numerical calculations of long-range dipo- 
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FIG. 1: (Color online) Using Gaussian D(Hs) with width 
a = 1.0 to calculate the M(H) curves and AH(M, AM) 
curves on a 2D square lattice with N = 100 2 hysterons and 
J cx = 0.4. (Top) J dp = 0. (Middle) J dp = 0.4. (Bottom) 
Jd P ~ 0.8. (Left) M(H) curves: main loop and 5 recoil curves. 
(Right) AH(M, AM) curves for the 5 recoil curves: (solid 
lines) numerical result; (dotted lines) mean-field approxima- 
tion. Here M (or AM) is normalized to the saturation value 
Ms = N and H (or AH) is normalized to the coercive field 
H c . 



lar interactions in a system with periodic boundary con- 
ditions, we utilized the efficient formalism described by 
Lekner [HI- 

For our numerical study of the AH(M, AM) method's 
reliability, we assume a Gaussian distribution D(H§) of 
width cr = l for a two-dimensional square lattice compris- 
ing of total N hysterons. Different system sizes ranging 
from 50 2 to 400 2 have been studied to estimate finite- 
size inaccuracies. Results presented here were calculated 
for 100 2 , which we found to be sufficiently precise in all 
cases. In our simulations, we normalize all parameters 
to a and vary both J cx and Jd p . For each parameter set 
(J cx , Jdp), we calculate the complete set of M(i?)-curves 
(both the saturation hysteresis loop and recoil curves), 
from which AH(M, AM) data sets are then extracted. 

The results displayed in Fig. 1 show several specific ex- 
amples for parameter sets (J cx , Jd P ) = (0.4,0), (0.4,0.4) 
and (0.4,0.8) . The simulated M(H) curves are shown 
in the left column. It is clearly seen that increasing 
the strength of dipolar interactions shears the hystere- 
sis loops substantially as expected. The right column 
displays the corresponding AH(M, AM) curves. The 
solid lines are the numerically extracted results from the 
simulated M(H) curves while the dotted lines denote 
the mean-field behavior according to the expression of 
AH G (M, AM). Comparing Fig. 1(d), (e) and (f), we find 
that dipolar interactions of intermediate strength make 
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FIG. 2: (Color online) Contour plots of the fit quality mea- 
sures as functions of J cx and Jd p . (a) Pd- (b) R 2 . Dotted 
lines indicate the range of Jd p for realistic materials. 



the system most mean-field like. 

To study the effect of dipolar interactions on the 
AH(M, AM) method's reliability in a more system- 
atic and quantitative way, we need to introduce quan- 
titative reliability measures. Obviously, the reliability 
range of the mean-field approximation, upon which the 
AH(M, AM) method is based, can be checked by means 
of a least-squares fit of AHq (M, AM) to the numeri- 
cal data and a subsequent analysis of the conventional 
fit-quality measures, such as: (1) the square of the mul- 
tiple correlation coefficient R 2 and (2) the percentage 
difference Pd between the fitting result and the input 
parameter a. By definition, R 2 = 1 and Pd = 
would correspond to perfect data fitting, i.e. the ex- 
actness of the mean- field limit. We therefore calculated 
a least-squares fit to AHq (M, AM) for the numerically 
extracted AH(M, AM) data for each parameter set of 
(J cx , Jdp)- From these fits, we then computed both Pd 
and R 2 . The results are shown in contour plots (see 
Fig.2). 

The shape of the contour plots is rather interesting. 
First, it is nearly symmetric along the diagonal direc- 
tion, i.e. Jdp/ J ox = 1- This clearly demonstrates that 
the roles of exchange and dipolar interactions in deter- 
mining D(Hq) are almost equally important. Individ- 
ually increasing either one will make the AH(M, AM) 
method less reliable, while increasing both of them with 
proper strength ratio of order 1 will substantially extend 
the reliability range. Second, the shape is not really sym- 
metric. It is tilted upwards and smoother on the high 
Jdp side than the high J ox side. This suggests that the 
AH(M, AM) method can clearly cope with higher dipo- 
lar interactions than exchange interactions, in agreement 
with previous micromagnetic tests @. 

The overall shape of the contours can be qualitatively 
explained by the interaction compensation effect. As we 
know, the intergranular exchange interactions are ferro- 
magnetic (FM) and short-range while the dipolar interac- 
tions are anti-ferromagnetic (AFM) and long-range. The 
competition between the two "opposite" interaction ten- 
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FIG. 3: (Color online) Correlations between the reliability 
measures: (a) Pd and r; (b) R 2 and r. Each data point 
represents a different parameter set (J ex , Jd P )- Data points 
with the same Jd P values are grouped and shown with the 
same point symbol, as indicated by the legend. 

dencies will yield a variety of system behaviors. Gen- 
erally speaking, as we steadily increase J ox from to 
higher values while keeping Jd p constant, we shift the sys- 
tem from the AFM-interaction-dominated regime to the 
mean-field regime and to the FM-interaction-dominated 
regime. Only within the interaction compensation re- 
gion, the FM and AFM interaction tendencies nearly can- 
cel each other. Consequently, the system is most mean- 
field like and the AH(M, AM) method becomes most 
reliable there. From the model Hamiltonian, we notice 
that if J ex / Jdp = 1, the exchange and dipolar interac- 
tions will cancel exactly for the nearest-neighboring hys- 
terons in the lattice. However, due to the long-range and 
distance-dependent features of dipolar interactions, this 
cancellation will not be exact for hysterons with longer 
distances. This explains why the interaction compensa- 
tion region is only roughly symmetric along the diagonal 
direction. 

To quantify the interaction compensation region or 
equivalently the reliability range of the AH(M, AM) 
method, one can define a critical value for each relia- 
bility measure, above which this method is sufficiently 
accurate. For example, we might define R 2 — 0.98 to be 
the critical value for R 2 in accordance with the best avail- 
able experimental data [gj. Then the reliability range 
for the parameter set ( J cx , Jdp) can be clearly seen from 
the region enclosed by the second highest contour in the 
R 2 plot. This particular contour will be referred to as 
the critical contour, inside which the system is virtu- 
ally mean-field like and the AH(M, AM) method is reli- 
able. Note that for practical recording media, the ratio of 
Jdp/c will probably be limited within the range of 0.2 ~ I 
as illustrated in Fig. 2. One can easily see that within the 
realistic Jd P /c range, the dipolar interactions improve 
the reliability range of the AH(M, AM) method up to 



higher J cx values. We also notice that for J cx /er < 0.2 
(read off from the intercept of the R 2 critical contour on 
the J ex -axis), higher dipolar interaction will only make 
the AH(M, AM) method worse. But this part is very 
small compared to the part where dipolar interactions 
make the AH(M, AM) method robust. 

Besides the fit-quality measures R 2 and Pd, there is 
a self-consistency-check measure, which is based upon 
data redundancy in between multiple recoil curves [?|. 
One can test data for deviations from this redundancy by 

means of a quantity r = ^ ^ ■ (rfj (M)) 2 where (M) 
is by definition identical to zero within the mean-field ap- 
proximation, so is r [l2j]. The specific advantage of this 
quantity r is that it can be directly calculated from data 
sets alone without the need for any data fitting. There- 
fore, it is important to analyze the possible correlation 
between the fit-quality measure (either R 2 or P^) and 
the deviation-from-redundancy measure r. Knowledge 
of this correlation will enable us to estimate the suit- 
ability of the AH(M, AM) method without any data fit- 
ting. Considering this, we calculated r from the numer- 
ical AH(M, AM) data for each parameter set of ( J ex , 
Jdp). Overall, the contour plot of r shows very similar 
features as of the ones being displayed in Fig. 2 for R 2 and 
Pd • To visualize and quantify the correlation between R 2 
(Pd) and r, we plot R 2 (Pd) vs. r for the complete set 
of different (J cx , J<j p ) parameter (see Fig. 3). We find 
that the data collapse fairly well onto a single line in the 
high R 2 or low \P^\ range, in which the utilization of the 
AH(M, AM) method is sensible and accurate. This in- 
dicates that R 2 and Pd are highly correlated with r in 
the regime where these quantities matter. Due to the 
knowledge of these correlations, one now has a criterion 
that enables a judgment on the usefulness and reliability 
of any AH(M, AM)-data set evaluation. For that, one 
simply determines the r value from experimental or mod- 
eling data sets, looks up the expected precision with the 
help of Fig. 3 and then decides if a further data analysis 
is warranted or not. 

In summary, we find that the presence of dipolar in- 
teractions similar in size to those of real perpendicular 
recording media makes the AH(M, AM) method sub- 
stantially more precise and robust. The deviation-from- 
redundancy measure r, which is a self-consistency-check, 
is found to be a good predictor of the AH(M, AM) 
method's reliability and can be utilized as a criterion to 
decide if a full scale data-analysis is warranted. 
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